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1. Introduction 

Experiments in past (SPS), present (RHIC) and future (LHC) colliders are attempting to recre- 
ate an early condition of the universe known as the quark-gluon plasma, where quarks and gluons 
exist without being bound into hadrons. Colliders explore the zero chemical potential region of the 
QCD phase diagram where lattice simulations indicate that a significant increase in the degrees of 
freedom happens above a certain critical temperature T c « 175 MeV (for a recent review see [jl]]). 

Heavy quarkonium dissociation has been proposed long time ago as a clear probe of the quark- 
gluon plasma formation in colliders through the measurement of the dilepton decay -rate signal [§]. 
Since higher excited quarkonium states are more weakly bound than lower ones, the expectation 
is that, as the temperature increases, quarkonium will dissociate subsequently from the higher to 
the lower states providing also a dynamical probe of the quark-gluon plasma formation (for some 
recent experimental data see 

In order to study quarkonium properties in a thermal bath at a temperature T, the quantity to 
be determined is the quarkonium potential V, which dictates, through the Schrodinger equation 



the real-time evolution of the wave function <I> of a QQ pair in the medium. In the full theory, 
V must come from a systematic expansion in l/m (non-relativistic expansion), the leading term 
being the static potential, and in the energy E (ultrasoft expansion). The potential will encode all 
contributions from scales larger than E and smaller than m. If the temperature lies in this range, 
the potential will depend on it, if the temperature is smaller than or of the same order as E, the 
potential will be temperature independent. 

The expansions in l/m and E are best implemented in QCD by means of effective field theories 
(EFTs), very much in the same way as this has been done in order to describe quarkonium physics at 
zero temperature [Qj. In the EFTs, the full dynamics will be more complicated than the Schrodinger 
equation ( |1.1[ ), since the EFTs will account both for potential and/or non-potential interactions. 
However, Eq. ( |1.1[ ) will provide the correct leading-order dynamics. 

In the last two years, there has been a remarkable progress in constructing EFTs for quarko- 
nium at finite temperature and in rigorously defining the quarkonium potential. In |6|], the static 
potential was calculated in the regime T S> 1/r > mo, where mo is the Debye mass and r the 
quark-antiquark distance, by performing an analytical continuation of the Euclidean Wilson loop 
to real time. The calculation was done in the weak-coupling resummed perturbation theory. The 
imaginary part of the gluon self energy gives an imaginary part to the static potential and hence a 
thermal width to the quark-antiquark bound state. In the same framework, the dilepton production 
rate for charmonium and bottomonium was calculated in [0, ||[|. In static particles in real-time 
formalism were considered and the potential for distances l/r~mo was derived for a hot QED 
plasma. The real part of the static potential was found to agree with the singlet free energy and the 
damping factor with the one found in [Qj. In [[It]], a study of bound states in a hot QED plasma was 
performed in a non-relativistic EFT framework. In particular, the hydrogen atom was studied for 
temperatures ranging from T <C ma 2 to T ~ m, where the imaginary part of the potential becomes 
larger than the real part and the hydrogen ceases to exist. An EFT framework in real time and weak 
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Figure 1: Quarkonium at finite temperature: energy scales and EFTs. 



coupling for quarkonium at finite temperature was developed in Qll| ] ; in the rest of the presentation, 
we will follow closely that approach. 



2. Scales and effective field theories 



Quarkonium in a medium is characterized by different energy and momentum scales; there are 
the scales of the non-relativistic bound state (v is the relative heavy-quark velocity): m, the heavy 
quark mass, mv, the scale of the typical inverse distance between the heavy quark and antiquark, 
mv , the scale of the typical binding energy or potential and lower energy scale, and there are the 
thermodynamical scales: the temperature T, the inverse of the screening length of the chromoelec- 
tric interactions, i.e. the Debye mass mo and lower scales, which we will neglect in the following. 

If these scales are hierarchically ordered, then we may expand physical observables in the ratio 
of the scales. If we separate explicitly the contributions from the different scales at the Lagrangian 
level this amounts to substituting QCD with a hierarchy of EFTs, which are equivalent to QCD 
order by order in the expansion parameters. At zero temperature the EFTs that follow from QCD 
by integrating out the scales m and mv are called respectively Non-relativistic QCD (NRQCD) 
and potential NRQCD (pNRQCD), see [||] for a review. We assume that the temperature is high 
enough that T ^> gT ~ mo holds but also that it is low enough for T <^m and l/r ~ mv > mp to 
be satisfied, because for higher temperature the bound state ceases to exist. Under these conditions 
some possibilities are in order. If T is the next relevant scale after m, then integrating out T from 
NRQCD leads to an EFT that we may name NRQCDhtl, because it contains the hard thermal 
loop (HTL) Lagrangian Jl2[]. Subsequently integrating out the scale mv from NRQCDhtl leads 
to a thermal version of pNRQCD that we may call pNRQCDhtl- If the next relevant scale after 
m is mv, then integrating out mv from NRQCD leads to pNRQCD. If the temperature is larger 
than mv 2 , then the temperature may be integrated out from pNRQCD leading to a new version 
of pNRQCDhtl- The hierarchies of scales that lead to these different EFTs are schematically 
illustrated in Fig. [[]. Note that, as long as the temperature is smaller than the scale being integrated 
out, the matching leading to the EFT may be performed putting the temperature to zero. 

In the following we will also assume that v ~ a s , which is expected to be valid for tightly 
bound states: Y(1S), J/y, ... . 
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The mass m is the largest scale in the system. This allows to integrate out m first and organize 
the EFTs as expansions in l/m. The leading order in the l/m expansion corresponds to the static 
limit of NRQCD: 

I n f 

S£ = -jF^F'^ + ^qiilpqi + ^iDoY + X^iDoX, (2-1) 
4 (=1 

where y (%) is the field that annihilates (creates) the (anti)fermion; qi merif light (massless) quark 
fields. Only longitudinal gluons couple to static quarks. The relevant scales in static NRQCD are: 
1/r, V, ... T, nip, ... . 

Since we are interested in the real-time evolution of the heavy quark-antiquark pair, it is con- 
venient to modify the contour of the partition function in order to allow for real times, see, for 
instance, Jl3|]. In real time, the degrees of freedom double, modifying the propagators into 2 x 
2 matrices. Despite this, the advantages are that the way in which calculations are carried out is 
very close to the one for T = EFTs, moreover, in the static quark sector, the second degrees of 
freedom, labeled "2", decouple from the physical degrees of freedom, labeled "1". The technical 
reason for this is that the [Sq (p)]n component of a static quark propagator vanishes, hence the 
unphysical static quark fields "2" never enter in any physical amplitude, i.e. any amplitude that has 
the physical fields "1" as initial and final states. It is also very convenient to chose the Coulomb 
gauge: in Coulomb gauge, only transverse gluons carry a thermal part, but they do not couple to 
static quarks. Finally, also the static quark-antiquark potential has a 2 x 2 matrix structure, which 
reads 



(2.2) 



In the following, whenever we speak about the potential, we mean the physical one, i.e. the entry 
V in the above matrix. 




3. Static quark antiquark at T < V 

If the temperature is very low, T < V, then it does not affect the potential, which may be 
derived by integrating out the scale 1/r from (|2~T|). This leads to pNRQCD in the static limit, 
whose degrees of freedom are quark-antiquark states (color singlet S, color octet O), low energy 
gluons and light quarks. The Lagrangian is organized as an expansion in r. 

& = - 7 + £ ?i + Tr { S 1 (ido ~ V s ) S + 0* (/Do - V ) O} 

4 (=1 

+V A Tr{O t r^ES + S t r- (? EO} + yTr{O t r^EO + O t Or-gE} + --- . (3.1) 

At leading order in r, the singlet decouples from the octet and its equation of motion is (ido— 
V S )S = 0. We may identify V s and V with the singlet and octet potentials. They are Coulombic: 
V s (r) = —C F — and V (r) = — ^— at leading order in a s (N c = 3, C F = 4/3). 
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Figure 2: Single and double lines stand respectively for color singlet and color octet quark-antiquark prop- 
agators. The curly line stands for the chromoelectric correlator (E" (t)<j> a i,(t ,0)E h (0)} , where (j> a b is a Wilson 
line in the adjoint representation, circles with cross stand for chromoelectric dipole interactions. 



Thermal corrections do not affect the potential, but affect the static energy and the decay width 
through loop corrections. The leading correction is carried by the diagram shown in Fig. ||. The 
real part of the diagram gives the following thermal correction to the static energy 



8E = \ Nc C F ^rT 2 f{N c aJ(2rT)) , 
3 K 



(3.2) 



dx — 

e x 
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where f(z) 



the diagram gives the thermal width 
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(3.3) 



where «b (&°) = 1/ {e k °^ T — 1) is the Bose statistical factor. Corrections coming from the scale m/j 
are suppressed by powers of mjj/T. The width T originates from the fact that thermal fluctuations 
of the medium at short distances may destroy a color-singlet QQ into an octet plus gluons. This 
process is specific of QCD at finite T. We will call this process the singlet to octet break-up 
phenomenon; in QCD the corresponding diagrams are shown in Fig. |3[ 




Figure 3: QCD diagrams responsible for the singlet to octet transition width in a thermal bath. 



In the limiting case T <C V, we have 
and r is exponentially suppressed. 



i;r^r 3 <E*(0).E*(0)) r . 



(3.4) 



4. Static quark antiquark atl/r>r>V 

In the situation 1/r S> T 3> V, integrating out T from pNRQCD modifies pNRQCD into a 
new EFT, pNRQCD H TL- With respect to pNRQCD, the Yang-Mills sector of the pNRQCD H TL 
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Lagrangian gets an additional hard thermal loop part [12], which modifies, for instance, the longi- 
tudinal gluon propagator at k° = as 



i 



T 



1117 



k 2 + m 2 D |k| {k 2 + m 2 D y 



(4.1) 





Figure 4: Diagrams contributing to the real part of the potential. The cross in diagram (a) means that we 
consider only one octet potential insertion into the free octet propagator; the shaded circle in diagram (b) 
stands for the gluon self-energy diagrams. 

Also the potential in pNRQCDhtl gets an additional thermal correction SV to the Coulomb 
potential of pNRQCD. The leading contribution to the real part of the color-singlet potential comes 
from the diagrams shown in Fig. |] and reads 

2 k 3 



Re 8VAr) - ^N C C F a 2 rT 2 - ^{3)C F ^ r 2 Tm 2 D + ^(3)N C C F a 2 r 2 T 3 . (4.2) 



The first term stems from diagram (a) and is of order g 2 r 2 T 3 x V/T, the other ones stem from 
diagram (b) and are of order g 2 r 2 T 3 x (mn/T) . 





Figure 5: Diagrams contributing to the imaginary part of the potential. The two crosses in diagram (a) mean 
that we consider two octet potential insertions into the free octet propagator. 

The leading contribution to the imaginary part of the color-singlet potential comes from the 
diagrams shown in Fig. || and reads 

Jm8V s (r) = -^L a 3 T +^ay-Tm 2 D (- + y E + -hw + \ -41n2 -2^) 



C(2)7 



+^-ln2N c C F a 2 r 2 T 3 . 



(4.3) 



The first term stems from diagram (a) and is of order g 2 r 2 T 3 x (V/T) 2 , the other ones stem from 
diagram (b) and are of order g 2 r 2 T 3 x (mo/T) 2 . The imaginary part of the diagram (a) may be 
traced back to the singlet to octet break-up phenomenon introduced above while the imaginary part 
of the diagram (b) is due to the imaginary part of the gluon self energy. This may be interpreted 
as due to the scattering of soft space-like gluons emitted by the heavy quarks with hard particles 
(gluons and light quarks) in the medium. In plasma physics, this phenomenon is known as Landau 
damping 
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Divergences appear in the imaginary part of the potential at order g 2 r 2 T 3 x (mo/T) 2 , which 
have been regularized in dimensional regularization (e = (4 — d)/2, where d is the number of di- 
mensions). They cancel in physical observables against loop corrections from lower energy scales. 
In order to illustrate the cancellation mechanism, let's consider the case l/r>T> mrj ^> V. 
Under this condition also the scale mo contributes to the potential. Integrating out mrj from 
PNRQCDhtl leads to an extra contribution 8V S to the potential coming from the diagram shown 
in Fig. U when the momentum flowing in the loop is of order mo and consequently the gluon prop- 
agator is taken to be the HTL resummed gluon propagator as in Eq. (|Q|). This extra contribution 
reads 

ReSV^-gWx^) 3 , (4.4) 
lm5V s (r) = -^La s r 2 Tm 2 D (--y E + \nK + \n^ T + 5 ^]. (4.5) 



The divergence in the imaginary part exactly cancels the one in (4.3). 

Summing up the real and imaginary parts of the potential corrections obtained from the scales 
T and mo, we end up with the thermal correction to the static energy 8E and the thermal decay 
width r respectively: 

k 3 oc 2 

8E = Nc C F a 2 rT 2 --i;(3)C F -r 2 Tm 2 D + -t;(3)N c C F a 2 r 2 T 3 , (4.6) 
9 2 K 3 



-2Im5V, = ^^a s 3 r 
3 s 



Cf 9 9 / T 2 C'(2)\ 871 , , 

-^a s r 2 Tm 2 D [2y E -fa— -l- 41n2 - 2^ J - — In 2 N C C F a 2 r 2 T 



(4.7) 



The (leading) non-thermal part of the static energy is the Coulomb potential —Cp(X s /r. The thermal 
width has two origins. The first term comes from the thermal break up of a quark-antiquark color 
singlet state into a color octet state. The other terms come from imaginary contributions to the 
gluon self energy that may be traced back to the Landau-damping phenomenon. The first one is 
specific of QCD, the second one would also show up in QED. Having assumed mp ^> V, the term 
due to the singlet to octet break up is parametrically suppressed by (V /mo) 2 with respect to the 
imaginary gluon self-energy contributions. The lnr 2 //n|, term is a remnant of the cancellation 
occurred between an infrared divergence at the scale T and an ultraviolet divergence at the scale 
m D . 



5. Static quark antiquark at T > 1 jr S> mo 

In the situation r>l/r> mp, integrating out T from static QCD leads to static NRQCDhtl> 
which, at one loop, is static NRQCD with the Yang-Mills Lagrangian supplement by the HTL 
Lagrangian. Subsequently, integrating out 1/r leads to a specific version of pNRQCDhtl where 
the Coulomb potential gets corrections from HTL insertions. The leading real correction comes 
from the diagram shown in Fig. ^, which gives 

Re8V s (r) = -^-a.rml . (5.1) 
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m D insertion 



Figure 6: Leading real thermal correction to the potential. The black dot stands for the insertion of the real 
part of the HTL gluon self energy. 



This is a correction proportional to a s /r x (rni£>) 2 . 

The leading correction to the imaginary part of the potential comes from the diagrams shown 
in Fig. [7| which give 



Im SV s (r) 



N 2 C F 



a s 3 r + ^-a,r 2 Tm 2 D 



+ y£ + in7r + ln(rju) 2 -l 



(5.2) 



The first term comes from diagram (a) in Fig. ^. It is proportional to a s /r x (rV) 2 x (Tr) and 
its origin may be traced back to the singlet to octet break-up phenomenon. The other terms come 
from diagram (b) in Fig. |7|. They are proportional to a s /r x (rnip) 2 x (Tr) and their origin may be 
traced back to the Landau-damping phenomenon. 



(a) 



i 

njOOOOrJOOOOOOrJ 
i 



(b) 



-inm 2 D T insertion 



Figure 7: Leading contributions to the thermal decay width. Here, the black dot stands for the insertion of 
the imaginary part of the HTL gluon self energy. 

Divergences appear in the imaginary part of the potential at order a s /r x (rm D ) 2 x (Tr). They 
cancel in physical observables against loop corrections from lower energy scales. In order to illus- 
trate the cancellation mechanism, let's consider the case T ^> \/r^> 3> V, which is similar to 
the one discussed in the previous section. Integrating out m/> from pNRQCDhtl leads to an extra 
contribution 8V S to the potential coming from the diagram shown in Fig. || when the momentum 
flowing in the loop is of order mo and consequently the gluon propagator is taken to be the HTL 
resummed gluon propagator as in Eq. ( |4.1[ ). These extra contributions are the same as in Eqs. ( p~4| ) 



and (|4.5D. 

Summing up the real and imaginary parts of the potential corrections obtained from the scales 
1/r and mj), we end up with the thermal correction to the static energy 8E and the thermal decay 
width r respectively: 



SE 



Cf 2 

-— a s rm D , 
r= ^LalT + ^a,r 2 Tm 2 D (-2y E -Hrm D f 



(5.3) 
(5.4) 



8 



Effective field theories for heavy quarkonium at finite temperature 



Antonio Vairo 



The (leading) non-thermal part of the static energy is the Coulomb potential — C F a s /r. Again the 
thermal width has two origins. The first term comes from the thermal break up of a quark-antiquark 
color singlet state into a color octet state. The other terms come from imaginary contributions to the 
gluon self energy that may be traced back to the Landau-damping phenomenon. Having assumed 
mj) S> V, the term due to the singlet to octet break up is parametrically suppressed by (V /ttio) 2 with 
respect to the imaginary gluon self-energy contributions. The \n(rttio) 2 term is a remnant of the 
cancellation occurred between an infrared divergence at the scale l/r and an ultraviolet divergence 
at the scale mo- 
lt is in the situation T ^> l/r ^ mo S> V that quarkonium in the medium melts, if we assume 
that the melting condition is ^binding ~ F, where ^binding is the quarkonium binding energy. Using 
the above results, the condition gives g 2 /r ~ g 2 Tm D r 2 In 1/ (trior). For l/r ~ tng 2 and mo ~ gT, 
this leads to the melting temperature r me i t i ng ~ mg 4 / 3 (lnl/g) -1 ' 3 , where, assuming g < 0.5, we 



have neglected In In l/g with respect to In \/g [ |10| , |14fl . 



6. Static quark antiquark at T > l/r ~ mo 



In the situation T S> l/r ~ mo, integrating out T from static QCD leads to static NRQCDhtl- 
Subsequently, both the scales l/r and mo have to be integrated out at the same time; this implies 
using HTL resummed gluon propagators in the matching procedure that leads to a new specific 
version of pNRQCDhtl- 



(a) 



(b) 



HTL propagator 



Figure 8: Diagram (a) shows the leading mass self energy contribution and diagram (b) the leading potential 
contribution to the static energy. Dashed lines stand for longitudinal HTL resummed gluon propagators. 



The real part of the static energy is provided at leading order by the two diagrams shown in 
Fig- §: 



E = Re [28m + 8V s (r)] = -C F Okm D 



(6.1) 



which is of order a s mo- The result is in agreement with early results on 8m and 8V S [|l5|, [T6| ]. 
The thermal decay width is provided at leading order by the three diagrams shown in Fig. |9|: 



N*C F 



■a. 



3 r 



2C F a s T 



dx 



sm(morx) 



rm D Jo (x 2 + l) 2 



(6.2) 



The first term is due to the singlet to octet break-up mechanism and is of order a s mo X (Vr) 2 x 
T /mo the other ones, which were first derived in [Q], are of order a s mo x T /mo S> cc & mo, i.e. 
larger than the real part of the energy (we recall that the binding energy is already of the same order 
as the decay width at the lower temperatures discussed in the previous section). The imaginary part 



of 8m is minus twice the damping rate of an infinitely heavy fermion [17] 
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(a) 



(b) 



(c) 



/ 



\ 



HTL propagator 



/ 



\ 



Figure 9: Diagram (a) is the leading diagram contributing to the singlet to octet break-up mechanism, 
diagram (b) contributes to the heavy quark damping rate and diagram (c) encodes the Landau-damping 
phenomenon. Dashed lines stand for longitudinal HTL resummed gluon propagators. 

7. Conclusions 

In a framework that makes close contact with modern effective field theories of non-relativi- 
stic bound states at zero temperature, we have discussed the real-time evolution of a static quark- 
antiquark pair in a medium of gluons and light quarks at finite temperature under the special as- 
sumption of weak coupling both for the non-relativistic and the thermal dynamics. For temperatures 
T ranging from values smaller to larger than the inverse distance of the quark and the antiquark we 
have derived the potential, the energy and the thermal decay width. 

The derived potential, V s , is neither the color-singlet quark-antiquark free energy nor the in- 
ternal energy (whose practical definition, at variance with the T = case, is plagued by many 
difficulties; for a recent critical discussion we refer to [|^]). It has an imaginary part and may 
contain divergences that eventually cancel in physical observables. 

The derived potential describes the real-time evolution of a quarkonium state in a thermal 
medium. At leading order, the evolution is governed by a Schrodinger equation. In an EFT frame- 
work, the potential follows naturally from integrating out all contributions coming from modes 
with energy and momentum larger than the binding energy. For T < V the potential is simply 
the Coulomb potential. Thermal corrections affect the energy and induce a thermal width to the 
quarkonium state; these may be relevant to describe the in medium modifications of quarkonium 
at low temperatures. For T > V the potential gets thermal contributions, which are both real and 
imaginary. 

Two mechanisms contribute to the thermal decay width: the imaginary part of the gluon self 
energy induced by the Landau-damping phenomenon, and the quark-antiquark color singlet to color 
octet thermal break up. Parametrically, the first mechanism dominates for temperatures such that 
the Debye mass mj) is larger than the binding energy, while the latter dominates for temperatures 
such that mo is smaller than the binding energy. Finally, it has been argued that quarkonium 
dissociation may be a consequence of the appearance of a thermal decay width rather than being 
due to the color screening of the real part of the potential; this follows from the observation that 
the thermal decay width becomes as large as the binding energy at a temperature at which color 
screening may not yet have set in. 
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